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Tu, mou Aec u soda, Kino o6sedem, a mo, Kau ckbo3hak, 
npoHUKaem e me6ji, mo ZAasoAem, a nmo oSuhak... 

H. BpodcKuu 

“You, my forest and water! One swerves, while the other shall spout 
Through your body like draught; one declares, while the first has a doubt. ” 

J. Brodsky 

“Mein Vater, mein Vater, und hoerest du nicht, 
Was Erlenkoenig mir leise verspricht?” 

“Sei ruhig, bleib ruhig, mein Kind! 

In duerren Blaettem saeuselt der Wind. ” 

J. W. Goethe 

CRIES AND WHISPERS IN WIND-TREE FORESTS 

VINCENT DELECROIX AND ANTON ZORICH 

To the memory of Bill Thurston 
with admiration for his fantastic imagination. 


Abstract. We study billiard in the plane endowed with symmetric Z 2 -periodic 
obstacles of a right-angled polygonal shape. One of our mane interests is de¬ 
pendence of the diffusion rate of the billiard on the shape of the obstacle. We 
prove, in particular, that when the number of angles of a symmetric connected 
obstacle grows, the diffusion rate tends to zero, thus answering a question of 
J.-C. Yoccoz. 

Our results are based on computation of Lyapunov exponents of the Hodge 
bundle over hyperelliptic loci in the moduli spaces of quadratic differentials, 
which represents independent interest. In particular, we compute the exact 
value of the Lyapunov exponent for all elliptic loci of quadratic differentials 
with simple zeroes and poles. 


1. Introduction 

The classical wind-tree model corresponds to a billiard in the plane endowed with 
Z 2 -periodic obstacles of rectangular shape; the sides of the rectangles are aligned 
along the lattice, see Figure [T] 


B BB B HxH H 

□ n □ n nNn 

□ □ m ■ ■ □ 

Figure 1. Original wind-tree model. 


The wind-tree model (in a slightly different version) was introduced by P. Ehren- 
fest and T. Ehrenfest [Ehj about a century ago and studied, in particular, by J. 
Hardy and J. Weber [HaWcl . All these studies had physical motivations. 

Several advances were obtained recently using the powerful technology of devi¬ 
ation spectrum of measured foliations on surfaces and the underlying dynamics in 
the moduli space. For all parameters of the obstacle and for almost all directions 
the trajectories are known to be recurrent j AH] ; there are examples of divergent 
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trajectories constructed in jD]; the non-ergodicity is proved in (FUj . It was proved 
in |T7 hT] that the diffusion rate is it does not depend either on the concrete 
values of parameters of the obstacle or on almost any direction and almost any 
starting point, see Figure [2] 


Figure 2. The diffusion rate | does not depend on particular 
values of the parameters of the rectangular scatterer: it is the same 
for the plane with horizontal walls having tiny periodic holes, for 
narrow periodic corridors between “chocolate plates” and for any 
other periodic billiard as in Figure [l] 

In other words, the maximal deviation of the trajectory from the starting point 
during the time t has the order of £3 for large t in the following sense: 

logdiam(trajectory for time interval[0, t]) 2 

i-s-oo log t 3 

Thus, this behavior is quite different from the brownian motion, random walk in 
the plane, or billiards in the plane with periodic dispersing scatterers: for all of 
them the diffusion has the order y/t (and, thus, the diffusion rate is |). 

We address the natural question “what happens if we change the shape of the 
obstacle?”. We do not have ambition to solve this problem in the current paper in 
the most general setting. We just plant the wind-tree forest with several interesting 
families of obstacles and study the diffusion rate as the combinatorics of the obstacle 
inside the family becomes more complicated. We show, in particular, that if the 
obstacle is a connected symmetric right-angled polygon as on Figure 0 then the 
diffusion rate in the corresponding wind-tree model tends to zero as the number 
of corners of the obstacle grows; see Theorem [j] for more precise statement. This 
result gives an explicit affirmative answer to a question addressed by J.-C. Yoccoz. 


Figure 3. The diffusion rate in this wind-tree forest tends to 
zero when the number of corners of the obstacle grows. 

Now, when we have showed that for certain species of wind-trees the sound in 
the wind-tree forest propagates “as a whisper” (in the sense that the diffusion rate 
tends to zero), the challenge is to prove that for certain other species it propagates 
“as a cry” with the diffusion rate approaching 1. 
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Question 1. Are there periodic wind-tree billiards with diffusion rate arbitrary 
close to 1? Are there continuous families of wind-tree billiards like this? What are 
the shapes of the obstacles which provide diffusion rate arbitrary close to 1? 


1.1. Strategy of the proof. We develop the approach originated in the pioneering 
work of V. Delecroix, P. Hubert, and S. Lclicvre [PHD who applied the results from 
dynamics in the moduli space to the wind-tree model. Several very deep recent 
advances in dynamics in moduli spaces are of crucial importance for us. 

Following [DHLI we reformulate the original billiard problem in terms of the de¬ 
viation spectrum for the leaves of directional measured foliations on the associated 
flat surface S. This part is quite elementary and straightforward. 

By recent deep result of J. Chaika and A. Eskin ICkEj (based on fundamental 
advances of A. Eskin, M. Mirzakhani, A. Mokhammadi |EMi| . |EMiMo| 1 almost 
all directions on every flat surface are Lyapunov-generic. Combining the tech¬ 
niques [Fb] , [Z] of deviation spectrum with the results from [DHL] we conclude that 
the diffusion rate of the wind-tree billiard as in Figure [3] equals to the Lyapunov ex¬ 
ponent A + {h*) = A + (i>*) of certain very specific integer cocycles h*,v* € iL 1 (S', Z), 
where the flat surface S is considered as a point of the orbit closure C(S) in the 
ambient stratum of meromorphic quadratic differentials; the vector space 7L 1 (5', C) 
containing h* and v* is considered as a fiber over S of the complex Hodge bundle 
over £(5), and the Lyapunov exponents are the Lyapunov exponents of the complex 
Hodge bundle over C(S) with respect to the Teichmiillcr geodesic flow. 

In the current paper we intentionally focus on the family of billiards as in Figure[3] 
since for this family the rest of the computation is particularly transparent. Namely, 
the flat surface S belongs to the hyperelliptic locus C := Q h yp[\ 2m ^ —l 2m ) over 
<2(l m , — l m+4 ). (In our particular situation, the hypere lliptic locus is, actually, 
“elliptic”: the genus of the covering surface is 1.) Moreover, applying the arguments 
analogous to those in lAl LZI , one immediately verifies that the family B of billiards 
is so large (in dimension) that it is transversal to the unstable foliation in the 
ambient invariant submanifold C , and, thus, for almost every billiard table n in B 
the orbit closure C{S) of the associated flat surface <S(n) coincides with the entire 
locus C. 

The flat surface S has genus one, so the complex Hodge bundle Hf~, = H_ j_ has 
single positive Lyapunov exponent A+. The fact that h,v are integer immediately 
implies that A + {h) = A + (u) = A) 1 ". 

Formula (2.3) from |EKZ2] expresses the sum Y^i=i °f positive Lyapunov 
exponents of H\_ in terms of the degrees of zeroes (and poles) in the ambient locus 
and in terms of the Siegel-Veech constant c area (£). Since in our particular case the 
genus of the surface is equal to one, we get a formula for the individual Lyapunov 
exponent Xf in which we are interested. 

Developing Lemma (1.1) from [ EKZ2] we relate the Siegel-Veech constant c area (C) 
of the hyperelliptic locus C := Q h yp(\ 2m y — l 2m ) over Q( l m , —l m=4 ) to the Siegel 
Veech constants cq(Q( l m , — l m + 4 )) of the underlying stratum in genus zero. Plug¬ 
ging in the resulting expression the explicit values of cc(Q( l m , —l m + 4 )) obtained in 
the recent paper |AtEZ| and proving certain combinatorial identity for the resulting 
hypergeometric sum we obtain the desired explicit value of A+(Q hyp (l 2m , — l 2m )), 
which represents the diffusion rate in almost every original billiard. 
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Remark 1. Our results provide certain evidence that when the genus is fixed and 
the number of simple poles grows, the Lyapunov exponents of the Hodge bundle 
tend to zero (see IGrHuj for the original conjecture). 

1.2. Structure of the paper. In section[5]we state the main result in two differ¬ 
ent forms. In section [3] we show how to reduce the problem of the diffusion rate in 
a generalized wind-tree billiard to the problem of evaluation of the top Lyapunov 
exponent of the complex Hodge bundle over an appropriate hyperelliptic locus 
of quadratic differentials. In section ld.ll we revisit the original paper rfTHLl where 
this question is treated in all details for the original wind-tree model with periodic 
rectangular scatterers. We suggest, however, several simplifications. Namely, in 
section [3721 we describe the hyperelliptic locus over certain stratum of meromorphic 
quadratic differentials in genus zero where lives the flat surface S corresponding to 
the wind-tree billiard and we show that the diffusion rate corresponds to the top 
Lyapunov exponent Af of the complex Hodge bundle over the PSL(2, R)-orbit clo¬ 
sure of S. Following an analogous statement in [AtEZl we prove in section [3731 that 
for almost any initial billiard table n the PSL(2, R)-orbit closure of the associated 
flat surface 5'(n) coincides with the entire hyperelliptic locus. At this stage we 
reduce the problem of evaluation of the diffusion rate for almost all billiard table 
in the family to evaluation of the single positive Lyapunov exponent exponent of 
the complex Hodge bundle H+ = H^ over certain specific hyperelliptic locus. 

In section [4] we evaluate this Lyapunov exponent. We start by recalling in sec¬ 
tion HU] the technique from IEKZ2| : we also relate the Siegel- Veech constant of the 
hyperelliptic locus with the Siegel-Veech constant of the corresponding stratum in 
genus 0. In section l4~2l we summarize the necessary material on cyclinder configura¬ 
tions in genus zero and on the related Siegel-Veech constants from [Bo] and [AtEZl . 
Finally, in section l4~3l we prove the key Theorem [2] evaluating the desired Lyapunov 
exponent Aj 1- . The proof uses a combinatorial identity for certain hypergeometric 
sum; this identity is proved separately in section [5] 

Following the title “ What’s next?'" of the conference, we discuss in appendix iBl 
directions of further research in the area relevant to the context of this paper. 


2. Main results 


Denote by J3(m) the family of billiards such that the obstacle has 4m corners with 
the angle n/2. Say, all billiards from the original wind-tree family as in Figures Q] 
and[2]live in B( 1); the billiard in Figure [3] belongs to S(3); the billiard in Figured 
belongs to H(17). 


Theorem 1. For almost all billiard tables in the family B(m) and for almost all 
directions the diffusion rate 5(m) is the same and equals 


S(m) 


(2m)!! 
(2m + 1)!! ' 


When m +00 d(m) has asymptotics 


8{m) 




Here the double factorial means the product of all even (correspondingly odd) 
natural numbers from 2 to 2m (correspondingly from 1 to 2m +1). For the original 
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wind-tree, when the obstacle is a rectangle, we have m = 1 and we get the value 
<5(1) = | found in [DHL] . 



Figure 4. The diffusion rate depends only on the number of 
corners of the obstacle and not on the particular values of (almost 
all) length parameters nor on the particular shape of the obstacle. 


Following the strategy described in section 11.11 we derive Theorem Q] from the 
following result. 


Theorem 2. The locus Q h vp (l 2m ,—l 2mS ) over Q(l m , — l m+4 ) is connected and in¬ 
variant under the action o/PGL(2,R). The measure induced on Q^ vp ( l 2m ,—l 2m ) 
from the Masur-Veech measure on Qi(l m , —l m+4 ) is PSL(2,R )-invariant and er- 
godic under the action of the Teichmuller geodesic flow. The Lyapunov exponent 
X+ (m) of the Hodge bundle H+ over the locus Q h vP( l 2m , — l 2m ) under consideration 
has the following value: 


( 2 . 1 ) 


A+(m) 


(2m)!! 
(2m + 1)!! 


Remark 2. Note that there is a very important difference between the case of m = 1 
(corresponding to the classical wind-tree) and the cases m > 2. Namely, the locus 
Q hyp (l 2 ,—1 2 ) over <2(1,— l 5 ) is nonvarying: the Lyapunov exponent A/~ = | for 
all flat surfaces in the locus Q h yp(l 2 , —l 2 ). For m > 2 it is not true anymore. First 
of all, for each integer m > 2, taking appropriate unramified covering of degree m 
of a flat surface in the stratum <2(1 2 , — l 2 ) we get a flat surface in the hyperelliptic 
locus of Q h vP( l 2m , —l 2m ) over <2(l m , —l m+4 ). Concretely, such surface can be built 
starting from the original wind-tree model with rectangles and taking a fundamental 
domain that is made of m copies of the unit square. 

By construction, the Lyapunov exponent A/~ of the resulting Teichmiiller curve 
in Q h yp( l 2m ,—l 2m ), and, hence, the diffusion rate 6 for the corresponding wind- 
tree billiard does not change: 5 = Xf = 2/3. Furthermore, for m = 2 we were able 
to find examples of square tiled surfaces for which the value is neither the generic 
value (5(2) = 8/15 = 0.5333 ... nor 2/3 = 0.6666 ... 


permutations r and u 


r = (1,2,3,4,5, 6 , 7)(8,9,10,11,12,13,14) 
u = (1,3,13,8,2,14) (4,6,11,5,10,12)(7,9) 

§ = 0.6060... 

r = (1,2,3,4,5, 6 , 7,8)(9,10,11,12,13,14) 
u = (1,2,3,14,9)(4,13)(5, 6,7,11,12)(8,10) 

£ = 0.5454... 


See also the example in Appendix [Al with m = 3. 
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Question 2. What are the extremal values of over all closed PSL(2, R)-invariant 
suborbifolds in a given stratum (given locus)? Same question for the concrete hy- 
perelliptic locus Q h yp(l 2m , — l 2m ) over Q(l m , — l m + 4 ) ? What are the shapes of 
billiards for which these values are achieved (if there are any wind-tree billiards 
corresponding to these invariant suborbifolds)? 

3. From billiards to flat surfaces 

3.1. Original wind-tree revisited. Recall that in the classical case of a billiard 
in a rectangle we can glue a flat torus out of four copies of the billiard table and 
unwind billiard trajectories to flat geodesics on the resulting flat torus. 

In the case of the wind-tree model we also start from gluing a flat surface out of 
four copies of the billiard table. The resulting surface is ZffiZ-periodic with respect 
to translations by vectors of the original lattice. We pass to the quotient over ZffiZ 
to get a compact flat surface without boundary. For the case of the original wind- 
tree billiard the resulting flat surface X is represented at Figure [5] It has genus 5; 
it belongs to the stratum TL{ 2 4 ) (see section 3 of [DHL for details). 

One of the key statements of [ DHL] can be stated as follows. 

Let n be the original rectangular obstacle, define the corresponding wind-tree 
billiard by the same symbol n. Let X = X(n) be the flat surface as in Figured] 
constructed by the wind-tree billiard defined by the obstacle n. 

Consider the SL(2, R)-orbit closure C(X) C TL( 2 4 ) of the flat surface X. Consider 
the cohomology classes h*,v* £ H 1 (X,Z) Poincare-dual to cycles 

h = hoo — hoi + h io — hn 
v = v 0 o - v w + v 01 - Vii 

(see Figure [5J as elements of the fiber over the point X £ £{X) of the complex 
Hodge bundle H over C(X). 

Theorem l IDHLj h The diffusion rate in the original wind-tree billiard n coincides 
with the Lyapunov exponent A (h*) = X(v*)) of the complex Hodge bundle Hwith 
respect to the Teichmiiller geodesic flow on C(X). 



Figure 5. The flat surface X obtained as a quotient over Z © Z 
of an unfolded wind-tree billiard table. 


Remark 3. Recent result in IC'kFI proving that for any flat surface almost all 
directions on it are Lyapunov-generic allows to simplify part of the argument in 
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the proof of the above Theorem. In particular, it justifies that A (h*) and X(v*) are 
well-defined for X endowed with almost all direction. 

Note that any resulting flat surface X as in Figure [5] has (at least) the group 
(Z/2Z) 3 as a group of isometries. As three generators we can choose the isometries 
t/j and t v interchanging the pairs of flat tori with holes in the same rows (corre¬ 
spondingly columns) by parallel translations and the isometry i acting on each of 
the four tori with holes as the central symmetry with the center in the center of 
the hole. 




Figure 6. A surface S in the hyperelliptic locus Q hyp (l 2 , — l 2 ) 

(on the left) is a double cover over the underlying surface S in 
<2(1, —l 5 ) (on the right) ramified at four simple poles represented 
by bold dots. 

Consider the quotient S of the flat surface S over the subgroup (Z/2Z) 2 of 
isometries spanned by and i o t v . The resulting surface S as on the left side 
of the Figure [G] belongs to the stratum Q(l 2 ,—l 2 ); in particular, it has genus 1. 
The surface S obtained as the quotient of the original flat surface X over the 
entire group (Z/2Z) 3 as on the right side of the Figure [S] belongs to the stratum 
Q(l, — l 5 ); in particular, it has genus 0. Clearly, S is a ramified double cover over 
S with ramification points at four (out of five) simple poles of the flat surface S. 

Lemma 3.1. Consider the natural projection p : X -A S . The following inclusion 
is valid: 

h* ep*{H\S;Z)). 

Proof. Recall that 

h = hoo — hoi + hio — hu , 

where the cycles hij, i,j = 0,1 are indicated in Figure [5] Now note that 

Th(hoo)=hi 0 L°T v (h 00 ) = -h 0 i l o t v o T h (h 0 o) = -hu 

Thus, the cycle h is invariant under the action of the commutative group T spanned 
by Tfi and to t v . This implies that the Poincare-dual cocycle h* is also invariant 
under the action of T and, hence, is induced from some cocycle in cohomology of 

S = X/T. □ 

Corollary 1. Consider the orbit closure A4(S) of the flat surface S'(II) £ Q( l 2 , — l 2 ) 
The diffusion rate in the original wind-tree billiard II in vertical direction coincides 
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with the positive Lyapunov exponent Af of the complex Hodge bundle H ^ with re¬ 
spect to the Teichmuller geodesic flow on A4(S). 

Proof. From the Theorem of Delecroix Hubert-Lelievre cited above we know that 
the diffusion rate coincides with the Lyapunov exponent X(h*). By the previous 
Lemma h* =p*(a) where a € iL 1 (5;Z). It is immediate to see that A (h*) = A(ct), 
where A(ct) is already the Lyapunov exponent of the of the complex Hodge bundle 
H £ with respect to the Teichmuller geodesic flow on M(S). Since g(S) = 1 the 
corresponding cocycle has Lyapunov exponents ±A+. Since a is an integer covector, 
X(a) cannot be strictly negative. Hence X(a) = A+. □ 

3.2. Elliptic locus and diffusion in the generalized wind-tree billiard. The 

hyperelliptic locus Q hyp ( l 2m , — l 2m ) in Q(l 2m , — l 2m ) is obtained by the following 
construction. For any flat surface S in Q(l m , — l m + 4 ) consider all possible quadru¬ 
ples of unordered simple poles. For each quadruple construct a ramified double 
cover with four ramification points exactly at the chosen quadruple of points. By 
construction the induced flat surface belongs to the stratum <2(l 2m , — l 2m ) in genus 
1. Considering all possible flat surfaces S in Q(l m , —l m + 4 ) and all covers over all 
quadruples of simple poles we get the locus Q hyp (l 2m , — l 2m ) over Q(l m , —l m+4 ). 

It is immediate to see that when the obstacle is a rectangle as in the original wind- 
tree billiard the surface S constructed in section 13.11 belongs to the hyperelliptic 
locus Q hyp (l 2 ,—l 2 ) over the stratum <2(1,— l 5 ). Similarly, when the obstacle has 
4 to corners with the angle 7r/2, the analogous surface S (as on the left of Figure 0 
belongs to the hyperelliptic locus Q%p(i 2 ™ ; _]2m) over Q(l m ,-l m + 4 ). 




Figure 7. A surface S in the hyperelliptic locus 

Q ft,yp (l 2m , _l 2m ) is a double cover over the underlying surface S 
in Q(l m , — l m + 4 ) branched at the four simple poles represented by 
bold dots. 

The arguments of Dclecroix-Hubert-Lelievre |, DHL] extend to the more general 
case of obstacles n symmetric with respect to vertical and horizontal axes with 
arbitrary numbers of angles \, basically, line by line with an extra simplification due 
to results |CkE| mentioned in Remark [3] Applying exactly the same consideration 
as above we prove the following statement generalizing Corollary [L] to arbitrary 
m £ N. 

Let n be a connected obstacle having 4m corners with angle j and 4(m — 1) 
corners with angle aligned in such way that all its sides are vertical or horizontal 
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(see the white domain in the left picture in Figure 0. Suppose that II is symmetric 
under reflections over some vertical and some horizontal lines. 

Proposition 1. Consider the orbit closure A i(S) C Q h yp(\ 2m , — l 2m ) 0 f f{ a t 
surface S'(n) in the ambient hyperelliptic locus Q h vP(l 2m ^ — l 2m ). The diffusion rate 
in the wind-tree billiard with periodic obstacles II aligned with the lattice coincides 
with the positive Lyapunov exponent of the complex Hodge bundle H ^ over M(S) 
with respect to the Teichmiiller geodesic flow on At (S'). 

3.3. Orbit closures. Consider the original wind-tree billiard with rectangular ob¬ 
stacles. Such a billiard is described by five real parameters: by two lengths of the 
sides of the external rectangle defining the lattice; by two lengths of the sides of the 
inner rectangle represented by the obstacle and by the angle defining the direction 
of trajectory, see Figure EO Varying continuously these parameters we obtain a 
continuous family of billiards. 

Starting with a more general obstacle as in Figure 0 having 4m corners of angle jj, 
4(m — 1) corners of angle ^ and symmetric with respect to vertical and horizontal 
axes of symmetry we get an analogous continuous family of billiards which we 
denote by B(m). It is immediate to check that 

(3.1) dimR B(m) = 2m + 2 

Here the direction of the billiard flow is not considered as a parameter of the 
family B. For the original wind-tree billiard with rectangular obstacles this gives 
diniR (13(1) x § 1 /(Z/2Z) 2 ) = 4 + 1 = 5, as we have already seen. 

Proposition 2. For Lebesgue-almost every directional billiard in H £ B(m) the 
GL(2,R )-orbit closure o/5'(n) in l m + 4 ) coincides with the entire ambient 

stratum Q(l m ,—l m+4 ). 

Proof. The Proposition is a straightforward corollary of Lemma |3.21 below. □ 

The following Lemma is completely analogous to Proposition 3.2 in |AtEZ| . 
Lemma 3.2. Consider the canonical local embedding 

B(m) x (§V(Z/2Z) 2 ) -a Q(l m ,-l m+4 ). 

For almost all pairs (n,0) in B(m) x (S 1 /(Z/2Z) 2 ) the projection of the tangent 
space T* (B(m) x (S 1 /(Z/2Z) 2 )) to the unstable subspace of the Teichmiiller geodesic 
flow is a surjective map. 

Proof. Let us first prove the statement for m = 1 and then make the necessary 
adjustments for the most general case. Consider a broken line following the upper 
part of the polygon as on the right side of figure 0 starting at the point Pq and 
finishing at the corner P 4 . Turn the figure by the angle 9. Consider the associated 
flat surface S € Q(l, — l 5 ). Consider the canonical orienting double cover S £ H(2). 
The six Weierstrass points of the cover correspond to the corners of our polygonal 
pattern and to the two bold points on the horizontal axe of symmetry. Thus the 
vectors of the resulting broken line considered as complex numbers are exactly the 
basic half-periods of the holomorphic 1-form u> corresponding to a standard basis of 
cycles on the hyperelliptic surface ( S ) (see sections 3.1 and 3.3 in |AtEZj for more 
details). The first cohomology U 1 (5 I ;C) serve as local coordinates in Q(l,— l 5 ). 
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The components of the projection of the vector of periods to the space i? 1 (5';M) 
are of the form 

±2sin(</>)|PjP i+ i| or ± 2 cos(</>)|PiPi + i| where i = 1,..., 4. 

Thus, for (j> different from an integer multiple of 7 t/ 2 the composition map 
T,(B(m) x (S 1 /(Z/2Z) 2 )) H\S-R) 

is a surjective map. 

It is clear from the proof that to extend it from m = 1 to arbitrary m £ N 
it is sufficient to show that the real dimension of B{m) coincides with the com¬ 
plex dimension of Q(l m , — 1 771+4 ). Recalling (13.11) and the classical formula for the 
dimension of a stratum Q(di,... ,d n ) of quadratic differentials 

dime Q(di, ..., d n ) = 2g + n — 2 

we conclude that 

diniR B(m) = 2 m + 2 = dime Q(l m , — f m + 4 ) 

which completes the proof of the Lemma. □ 

Combining Propositions [T| and [5] with the result of Chaika-Eskin ICkEj telling 
that for any flat surface almost all directions are Lyapunov-generic, we obtain the 
following Corollary, which proves the first part of Theorem [T] 

Corollary 2. For Lebesgue-almost every directional billiard in II £ B(m) the dif¬ 
fusion rate in almost all directions 6 6 (§ 4 /(Z/2Z) 2 ) in the wind-tree billiard II 
coincides with the positive Lyapunov exponent A + of the complex Hodge bundle 
H £ over the hyperelliptic locus Q h yp(l 2m , — l 2m ) w nh respect to the Teichmuller 
geodesic flow on this locus. 

It remains to evaluate the Lyapunov exponent A + of the complex Hodge bundle 
H £ over the hyperelliptic locus Q hyp (l 2m , —l 2m ), which we do in the next section. 

4. SlEGEL-VEECH CONSTANTS AND SUM OF THE LYAPUNOV EXPONENTS OF THE 

Hodge bundle over hyperelliptic loci 

In this section we relate the Siegel -Veech constants of an invariant hyperelliptic 
locus and of the underlying stratum of meromorphic quadratic differentials with at 
most simple poles on CP 1 . Applying the technique from lEKZ2j and developing the 
results from [ AtEZj this allows us to get an explicit value for the desired Lyapunov 
exponent A + of the hyperelliptic locus and thus, to prove Theorem [5] 

4.1. Sum of the Lyapunov exponents of the complex Hodge bundle over 
hyperelliptic loci of quadratic differentials. We need the following result 

Theorem (Theorem 2 in 1 lbKZ21 i. Consider a stratum Qi(d \,..., d n ) in the moduli 
space of quadratic differentials with at most simple poles, where di+- ■ H -d n = 4^—4. 
Let Aii be any regular PSL(2, R)-invariant suborbifold of Qi(di ,... ,d n ). 

The Lyapunov exponents A + > • • • > A + of the complex Hodge bundle H+ = H £ 
over A4i along the Teichmuller flow satisfy the following relation: 


^ T 3 C a rea (AA 1 ) 


(4.1) 


A + -I-+ A + 
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where 


(4.2) 



and c area (Mi) is the Siegel-Veech constant corresponding to the suborbifold Mi. 
By convention the sum in the left-hand side of equation (ED is defined to be equal 
to zero for g = 0. 

In the context of this paper we are particularly interested in the case when M 
is a hyperelliptic locus over some stratum of meromorphic quadratic differentials 
with at most simple poles in genus zero. 

Let Mi be a closed SL(2, R)- (correspondingly PSL(2, R)-invariant) suborbifold 
in a stratum of Abelian (correspondingly quadratic) differentials; let v be the asso¬ 
ciated SL(2, R)-ergodic (correspondingly PSL(2, R)-ergodic) measure on M\. 

Consider a locus M i of all possible double covers of fixed profile over flat surfaces 
from Mi. Suppose that it is closed, connected and SL(2,R)- (correspondingly 
PSL(2,R)-invariant), and that v is the ergodic measure on M such that v is the 
direct image of v with respect to the natural projection Mi —> Mi. 

Let cc be an area Siegel-Veech constant associated to the counting of multiplicity 
one configuration C of cylinders weighted by the area of the cylinder. (Here the 
notion “configuration” is understood in the sense of |MsZ| and |W2j : “multiplicity 
one” means that C contains a single cylinder). The reader might think of M 
as of some stratum of quadratic differentials in genus 0; then there are only two 
types of configurations (see section IT^l below or the original paper jBoj ). and both 
configurations contain a single cylinder. 

We assume that the profile of the double cover does not admit branching points 
inside the cylinders of the configuration C. Then the configuration C induces a 
cylinder configuration C on the double cover. Let c be the associated area Siegel 
Veech constant. The Lemma below relates c and c. it is a slight generalization of 
Lemma 1.1 in |EKZ2] . 

Lemma 4.1. If the circumference of the cylinder in the configuration has nontrivial 
monodromy (that is if the lift of the cylinder in the double cover is a unique cylinder 
twice wider), then c = c/2. 

If the circumference of the cylinder has trivial monodromy (that is, if the preim¬ 
age of the cylinder consists in two cylinders isometric to the one on the base), then 


c = 2c. 


Proof. The second case corresponds to Lemma 1.1 in |EKZ2| . The first case is 


proved analogously. 


□ 


In the same setting, let k be the expression (14.21) in degrees of singularities of the 
flat surfaces in the invariant manifold M and R be analogous expression in degrees 
of singularities of the double covers of fixed profile in the invariant manifold M . 

It would be convenient to introduce the following notations 


Ak := k — 2k 


^CarecL ■ — Carea(Mll) 2c ar ea(Mt l) . 


Lemma 4.2. For any ramified double covering the degrees of zeroes of the quadratic 
differential on the underlying surface and of the induced quadratic differential on 
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the double cover satisfy the following relation: 


(4.3) 


Ak = k — 2ft 


- E • 

4 ^ dj + 2 

ramification J 
points 


Proof. The non ramified zeros cancel out. For the other zeroes (and poles) a sin¬ 
gularity of degree dj at the ramification point gives rise to a zero of degree 2 dj + 2. 
Hence, 


: —2k = — 

24 ^ 

ramification 

points 

= - E 


(2 dj + 2)(2 dj + 6) _ ^dj(dj + 4 


2 dj + 4 


dj H- 2 


ramification 

points 


(dj + + 3) — dj(dj + 4) 

dj 2 


- E 

io z^ 


12 'r—' ^ + 2 

ramification J 
points 


□ 


The following notational Lemma would help to simplify certain bulky computa¬ 
tions. 

Lemma 4.3. For any locus Mi of double coverings of fixed profile as above over 
a PSL(2, M)-invariant orbifold Mi in some stratum of meromorphic quadratic dif¬ 
ferentials with at most simple poles on CP 1 the sum of Lyapunov exponents 

A+ = Xf + ... A+ 

satisfies the following relation 

7T 2 

(4.4) A+ = Ak + — ■ A c area . 

Proof. For any invariant submanifold A4i in a stratum of meromorphic quadratic 
differentials with at most simple poles on CP 1 the sum of Lyapunov exponents is 
null, so formula (EH) gives 

7r 2 

0 = ^ T • c area (Mi). 

By the same formula (14.11) gives we have 

7T 2 

A — k T • c area (Mi). 

Extracting from the latter relation twice the previous one we obtain the desired 
relation (14.41) . □ 

4.2. Configurations for the strata in genus zero and corresponding Siegel 
Veech constants (after Boissy and Athreya Eskin Zorich). In this section 
we recall briefly the results from |Bo] describing configurations of periodic geodesics 
for flat surfaces in genus zero, and the results from |AtEZj providing the val¬ 
ues of the corresponding Siegel Veech constants. By Q(di,... ,dk) we denote a 
stratum of meromorphic quadratic differentials with at most simple poles, where 
di G {—1,1,2,...} denote all zeroes and poles, and Xa=i di = — 4. 
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A “pocket”. In this configuration we have a single cylinder filled with closed 
regular geodesics, such that the cylinder is bounded by a saddle connection joining 
a fixed pair of simple poles Pj 1 , Pj 2 on one side and by a separatrix loop emitted 
from a fixed zero Pi of order d,; > 1 on the other side. 

By convention, the affine holonomy associated to this configuration corresponds 
to the closed geodesic and not to the saddle connection joining the two simple poles. 
(Such a saddle connection is twice as short as the closed geodesic.) 



Figure 8. A “pocket” configuration with a cylinder bounded on 
one side by a saddle connection joining two simple poles, and by a 
saddle connection joining a zero to itself on the other side. 


By Theorem 4.5 and formula (4.28) in lAtEZI . the Siegel-Veech constant 
corresponding to this configuration has the form 


(4.5) 


vocket _ di T 

c h,h;i - ( fc _ 4) ' ^2 ■ 


One can consider the union of several configurations as above fixing the pair of 
simple poles Pj 1 , Pj 2 but considering any zero Pi on the boundary of the cylinder. 
By Corollary 4.7 and formula (4.36) in [AtEZj . the resulting Siegel-Veech constant 
c-t 4 corresponding to this configuration has the form 


(4.6) 


pocket 
'3i ,32 


l 

2-7T 2 


A “dumbbell”. For the second configuration we still have a single cylinder 
filled with closed regular geodesics. But this time the cylinder is bounded by a 
separatrix loop on each side. We assume that the separatrix loop bounding the 
cylinder on one side is emitted from a fixed zero Pi of order d,; > 1 and that the 
separatrix loop bounding the cylinder on the other side is emitted from a fixed zero 
Pj of order dj > 1. 

Such a cylinder separates the original surface S in two parts; let 
be the list of singularities (zeroes and simple poles) which get to the first part and 
Pj 1 ,..., Pj k2 be the list of singularities (zeroes and simple poles) which get to the 
second part. In particular, we have i £ [ii,... ,i^} and j £ {ji, ■ ■ ■ ,jk 2 }- We 
assume that S does not have any marked points. Denoting as usual by dk the order 
of the singularity Pk we can represent the sets with multiplicities a := {di,..., dk} 
as a disjoint union of the two subsets 

{di; ■ • • > dk} = {djj,... di ki } U {dj 1 ,..., dj k2 }. 

(Recall that {di,..., dk} denotes all zeroes and poles.) This information is consid¬ 
ered to be part of the configuration. 
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Figure 9. A “dumbbell” composed of two flat spheres joined by 
a cylinder. Each boundary component of the cylinder is a saddle 
connection joining a zero to itself. 


By Theorem 4.8 and equation (4.38) in (AtEZl the corresponding Siegel-Veech 
constant c dumbeU is expressed as follows: 

(A 7] dumbell _ ( d * + l )( d i + X ) (fci - 3)! (fc 2 - 3)! 1 

{ 2 ‘ (k — 4)! 'tt 2 ' 

According to [Bo] and [MsZj . almost any flat surface S in any stratum <2i(di,..., dk) 
of meromorphic quadratic differentials with at most simple poles different from the 
pillowcase stratum Qi(—l 4 ) does not have a single regular closed geodesic not con¬ 
tained in one of the two families described above. 

Finally, by Corollary 4.10 from |AtEZ| (generalizing the theorem of Vorobets 
from [Vo] ) the Siegel-Veech constant c area is expressed in terms of the above Siegel 
Veech constants as follows. For any stratum Qi(di,... ,dk) of meromorphic qua¬ 
dratic differentials with simple poles on CP 1 the Siegel-Veech constant c ar ea is 
expressed in terms of the Siegel-Veech constants of configurations as follows: 

1 \ - 

(4.8) Carea ~ ^ , CC ■ 

Configurations C 
containing a cylinder 

4.3. Proof of Theorem [2} Now everything is ready for the proof of Theorem [2] 

Proof. The connectedness of our hyperelliptic locus follows from the fact that it 
can be seen as a C*-bundle over the space C(m, m, 4) of configurations of points on 
CP 1 considered up to a modular transformation. More precisely, C(m, to, 4) can be 
seen as a set of configurations of 2m + 4 distinct points arranged into three groups 
(“colored into three colors”) of cardinalities {to, to, 4}, where the points inside each 
groups are named. The first group represents the simple zeroes, the second one 
unramified simple poles, and the last one — the ramified simple poles. Clearly, 
such space is connected. 

By the results of H. Masur |Ms| and of W. Veech IVej . the Teiclnniiller geodesic 
flow on the underlying stratum is ergodic, and, moreover, “sufficiently hyperbolic”. 
Thus, the induced flow on any connected finite cover is also ergodic. The hyperel- 
liptic locus Q^P(l 2m , _i 2m ) over <2(l m ,— l m + 4 ) is a finite cover, and as we have 
just proved it is connected. This proves the ergodicity of the Teichmiillcr flow on 
the hyperelliptic locus. 

Let us compute Ad™ c e k a et , see section ixn When both simple poles Pj ± , P n in¬ 
volved in the “pocket” configuration are unramified points or when they are both 
ramified, the holonomy of the hyperelliptic cover S —> S along the perimeter of 
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the cylinder is trivial, so by Lemma 14.11 such configurations do not contribute to 
Ac^°ea e * ■ By Lemma 14.11 a configuration when one of P 3l , Pj 2 is ramified and the 
other one is nonramified, contributes to A^ ei with a weight — |. Since we have 4 
ramified poles and to nonramified, there are 4m such configurations. Applying (14.61) 
and (HU) with k = m + (m + 4) we get 


(4.9) 


• Aff^ et = — • 4m ■ - 


1 


2 to + 1 


1 

2-7T 2 


m 


2 to + 1 


Let us proceed to computation of Ac„™ bben . When the number of ramified 
simple poles on two parts of the “dumbbell” is even, the holonomy of the cover 
along the perimeter of the cylinder is trivial, so by Lemma 14. 1 1 such configurations 
do not contribute to Acf“™ bbeii . By Lemma IXTI a configuration when the number of 
ramified simple poles on each side of the dumbbell is odd contributes to Ac)(“™ bben 
with a weight — To compute the number of such configurations we remark that 
we have to split m named simple zeroes into two groups of mi and to — TOi ones; 
we also have to split 4 ramified simple poles into 1 and 3; finally we have to split 
to unramified simple poles into mi + 1 and to — mi — 1 ones to have in total mi + 2 
simple poles on one side of the “dumbbell” and to — TOi + 2 on the other side. Note 
that the fact that there is a single ramified simple pole on one part and 3 ramified 
simple poles on the other makes our count of configurations asymmetric. Finally 
note that we have to chose one of m\ zeroes to be located at the boundary of the 
cylinder on one side and one of m — mi zeroes to be located at the boundary of the 
cylinder on the other side. 

For any given mi, where 1 < mi < to — 1 our count gives 



• (toi • (to — mi)) 


Applying the general formulae (1X71) and (14.81) to cf“™ be e “; taking into consid¬ 
eration that in our particular case we have di = dj = 1; k\ = m\ + (mi + 2); 
= (to — mi) + (m — mi + 2); and k = m + (to + 4); and applying (|4.8D we get 


(4.10) yAc“ J 

2 m —1 / 

= iE 


TO 

TOi 


TO 

TOi — 1 


(toi • (m — TOi)) 


(1 + 1)(1 + 1) ((2mi+2)-3)!-((2m-2mi + 2)-3)! 1 


(2m + 4) - 3 


2 TO 


E 




+ 1 \mij Vtoi — 1 

7Tli=l x X X 


m 


(2m+ 4)-4)! 

(2mi)! • (2 to — 2toi)! 


(2m)! 

771—1 / 771 W 771 \ 

^ ^ \771i / \7711 — 1/ 

771i=l 


2ra+ j ^ (|“) 


Summing up (14.91) and (14.101) and applying the standard convention 

:= 0 and 


TO 

to + 1 


:= 1 
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we obtain 

1 (mi) Lr+l) = _l< ( 2w ) !! 

2m+1 ^ ( 2m ) (2m+ 1)!!' 

mi=0 \2m 1 J v ' 

where the second equality ils the combinatorial identity (15.21) proved in Proposi¬ 
tion [3] 

By formula (14.31) we have 


(4.11) 


~ ^C ar f>a 


(4.12) 




4 £-(-1 + 2 ) 


= 1 


Plugging the results (14.111) and (14.121) of our calculation in the general for¬ 
mula (14.41) for the sum A + we get 


A + = 1 + 



(2m)!! \ 

(2m + 1 )!!) ' 


It remains to note that the stratum Q(l 2m ,—l 2m ) corresponds to genus one, so 
the spectrum of Lyapunov exponents of H+ contains a single entry and A + = A) 1 ", 
which proves Theorem [2] □ 


Corollary 3. The Lyapunov exponent (m) tends to zero as m tends to infinity. 
More precisely, 


< 4 ' 13 > ■ s <"'> = ^=( 1 + 0 (£>)' 

Proof. Rewriting the double factorials in terms of usual factorials as 

(2m)!! = (( 2 m)\\) 2 = 1 _ 2 2 m (m|f 

(2m+ 1)!! (2m+1)! 2m+1 (2m)!’ 

and applying the Stirling’s formula 



to both factorials and simplifying the resulting expression we get 


(2m)!! 
(2m + 1)!! 


2m + 1 



Clearly, the latter expression tends to zero as m tends to infinity. 


□ 
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5. Combinatorial identities 
Proposition 3. For any ro€ff the following identities hold 


f-fflg) (2")'' ,m ( m! ) 2 

ti CS (2™- 1 )'! (2m)l 


(5.2) 


m /m\ f m \ 

E \k) Vfc+1/ 
(2 m\ 
k =0 V2 k) 


2m + 1 — 


(2 to )!! 
(2m — 1)!! 


(5.3) 


E 


fc=o 



= 2m + 1. 


Proof of Proposition 0 First note that 


m + 1\ fm+ 1\ /m 

fc + iy = V y + V fc 


Thus, the expression in the left-hand side of (15.311 is the sum of the corresponding 
expressions in the left-hand sides of (15.11) and (15.21) . Hence, any two out of three 
identities (15.11) (15.31) imply the remaining one. 

Proof of identity (15.31) . Developing binomial coefficients into factorials in (15.31) 
and moving the common factorials to the right-hand side, we can rewrite this iden¬ 
tity as 


(5.4) 


E 


k—0 


(2fc)! 

*!(* + !)! 


(2m — 2fc)! 

(m — ft)! (m — ft)! 


cr)- 


Denote the sum in the left-hand side of (15.41) by S 3 (m). We prove the latter identity 
by induction. For m = 0 it clearly holds. Assuming that (15.41) holds for some integer 
m we are going to prove that 


(5.5) 


(m + 2) • S 3 (m + 1) — 2(2m + 3) • s 3 ( 171 ) = 0 . 
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Since the right-hand side of (15.411 satisfies the same relation (which is an exercise), 
this completes the step of induction. It remains to prove (15.511 under assump¬ 
tion (15.41) . 


(m + 2) • sz{m + 1) — 2(2 m + 3) • S 3 (m) = 


m+1 


= (m + 2) ^ 


( 2 * 0 ! 


(2m + 2 — 2fc)! 


k =0 


k\ {k + 1)! (m + 1 — k)\ (m + 1 — k)\ 


— 2(2m + 3) ^ 


(2 k)\ 


(2m — 2 k)\ 


fc=o 


fc! (fc + 1)! (m — k)\ (m — k)l 


— (m + 2) 


(2m+ 2)! 


£ 


(2 fc)! 


(2m — 2 fc)! 

(m + l)!(m + 2)! ' fc! (fc + 1)! (m — A;)! (m — fc)! 


, . (2m + 2 — 2fc)(2m + 1 — 2 k) , . , 

(m + 2)^——;- + | i -^ - 2(2m + 3) = 


= 2 ( 2 ^+l) ! 


(m + 1 — k)(m + 1 — k) 

(2 fc)! (2m — 2fc)! 


(m+l)!m! “+!(fc + l)! (m — fc)!(m — fc)! 


- 2 - 


fc + 1 


m + 1 — k 


= 2 


= 2 


2 m + 1 
m 

2m + 1 
m 


-2£ 


(2m- 2fc)! 


fc—o 


(2fc)! _ 

fc!fc! (m — fc)! (m — fc + 1)! 


-2£ 


3=0 


(2(m — j))! _ 

(m-j)!(rn-j) ! j!(j + l)! 


(2j)! 


= 2 


2m + 1 
m 


— 2sz{rn) = 0 . 


where the last equality is the induction assumption. Identity (15.41) and hence iden¬ 
tity (15.31) is proved. 

Proof of identity (15.11) . Developing binomial coefficients into factorials in (15.11) 
and simplifying common factorials in the right and left hand side of EU , we can 
rewrite this identity as 

'2j\ ( 2m — 2 ? n 

JU-j 


E 

3=0 


= r 


which is identity (3.90) in IGoj . 
Proposition [3] is proved. 


□ 


Appendix A. Removing some squares in the wind-tree model 

Let us consider the periodic wind-tree models with square obstacles. If we remove 
periodically one obstacle out of four, we can still perform the same construction as 
before and end up with a surface in Q h yp( l 6 , —l 6 ) over Q(l 3 , — l 7 ). 

Namely, unfolding the wind-tree in Figure [TO] and taking the quotient over Z®Z 
we get a compact translation surface represented in Figure 111! This translation 
surface corresponds to the wind-tree in Figure [TO] exactly in the same way as the 
compact flat surface in Figure [5] corresponds to the original wind-tree in Figure [l] 

As in the previous examples, the flat surface X in Figure [IT] has the group 
(Z/2Z) 3 as a group of isometries (we have already seen exactly the same group of 
symmetries for the surface in Figure [5]). As before, we can choose as generators 
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Figure 10. A square-tiled wind-tree where we regularly remove 
one obstacle in every repetitive pattern of four. 


F 




1 




1 

> 


i 


r 









T 




Figure 11. The translation surface X £ TL( 2 12 ) obtained as the 
quotient over Z © Z of an unfolded wind-tree billiard table in the 
left picture in Figure fT(Tl 

the isometries and t v interchanging the pairs of flat tori with holes in the same 
rows (correspondingly columns) by parallel translations and the isometry i acting 
on each of the four tori with holes as the central symmetry with the center in the 
center of the square. Passing to the quotient over (Z/2Z) spanned by t v we get the 
square-tiled surface S £ "H(2 6 ) as in Figure [T^l 
It is encoded by the following two permutations 

r = (1, 2, 16. 14.15, 3)(4, 5, 6, 7)(8, 22)(9, 21)(10, 11, 12,13)(14, 15,16)(17, 18, 19, 20)(23, 24, 25, 26) 
u = (1,4,10)(2, 5, 9,11)(3, 7, 8,13)(6, 12)(14, 17, 23)(15,18, 22, 24)(16, 20, 21, 26)(19, 25). 

We construct its SL(2, K)-orbit; analyze its elements, and apply formula (2.12) 
from IEKZ2I to find 

oacq 

(A.l) A, (£) + ■••+Aren¬ 

as the sum of the Lyapunov exponents of the resulting arithmetic Teichmiiller curve 
C(S). 
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Figure 12. The flat surface S = X/t v . 


Now, following the strategy of section [3~il we pass to the second quotient S = 
X/(rh,T v ) £ Q h vp( l 6 , —l 6 ) (on the left of Figure flfll) followed by the quotient 
S = X/{Th,T v ,i) £ Q(l 3 ,—1 3+4 ) (on the right of Figure IbH) . 



Figure 13. The flat surfaces S = X/{Th,r v ) £ Q hyp (l 6 , — l 6 ) 

(on the left) and S = X/{Th,r v ,i} £ Q(l 3 ,—1 3+4 ) (on the right). 

We inverse shadowing at this picture with respect to Figures IT2l 
and EH now we shadow the surface, and not the obstacles. Tiny 
black discs at the vertices represent the simple poles; the other 
vertices correspond to simple zeroes. 

The surface S £ T~L( 2 6 ) is the orienting double cover of the surface S. Hence, the 
Lyapunov spectrum (IA.1I) of the arithmetic Teichmiiller curve C(S) is the union of 
the spectra of Lyapunov exponents of the arithmetic Teichmiiller curve C(S): 

{A i(S), ..., A 7 (5)} = {Ar(S),..., Ae (S)} U {A +(S)} . 

By formula (2.4) from |EKZ2j we have 

( Ar (S) + ... + A 6 -(S))-A+ = r ^ 1 —^=l.6.(± + lj=2. 

j such that J 
dj is odd 

Together with (1 A. 1 1> this implies that 

1 f 3088 \ 491 

2 \ 1053 J ~ 1053' 

and, hence, that the diffusion rate for the square-tiled wind-tree as in in Figure flOl 
is given by 

491 
1053' 


A t(S) 


5 = \+(S) 
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We shall see that in all the other cases, removing some of the obstacles out of 
every repetitive block of 2 x 2 obstacles in the wind-tree model with any parameters 
a and b we keep the diffusion rate 6 = Xf = 2/3. There are three cases to consider: 

(1) removing two obstacles which are in the same row or in the same column; 

(2) removing two obstacles which are on the same diagonal; 

(3) removing three obstacles. 


In the first and in the last case we can just choose a new fundamental domain of 
the rectangular lattice (duplicating it in the first case and choosing it 2 x 2 bigger 
in the third case) to reduce the situation to the original wind-tree with different 
parameters. We have seen that the diffusion rate in the original wind-tree as in 
Figure Q] does not depend neither on the parameters of the lattice, nor on the 
parameters of the obstacle. Hence in both cases S = A/" = 2/3. 

Let us consider the case where we remove two obstacles on the same diagonal as 
in Figure HU In this case, we can modify the construction a little bit. We unfold 
the billiard as before, but then we quotient the resulting periodic surface by the 
integer sublattice spanned by the integer vectors (1,1) and (1,-1) in Z © Z and 
not buy the entier lattice Z © Z as before. The quotient belongs to the same locus 
as the original wind-tree. We deduce that we again have <5 = A/" = 2/3. 

Let us emphasize that the above construction is very specific to the sublattice 
Z(l, 1) © Z(l, —1) which is invariant under reflexions by the horizontal and vertical 
axes. In such situation the quotient keeps a (Z/2Z) 3 group of symmetry. 



■ 


Figure 14. A square-tiled wind-tree where we regularly remove 
every second square obstacle has the original diffusion rate S = 2/2,. 


Appendix B. What is next? 

The name of the conference and personality of Bill Thurston suggest to discuss 
what are the current directions of research in the area. We focus on those aspects 
which are somehow related to the current paper. The selection represents our 
particular taste and might be subjective. These problems are, certainly, known in 
our community; some of them are already subjects of intensive investigation, so we 
do not claim novelty in this discussion. 

The Magic Wand Theorem of Eskin Mirzakhani -Mohammadi enormously am¬ 
plified the importance of the classification of the orbit closures. Ideally, any problem 
about an individual flat surface (say, diffusion rate for a wind-tree model with ob¬ 
stacles of some specific “rational shape”) should be solved as follows: touch the 
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corresponding surface with the Magic Wand and find the corresponding orbit clo¬ 
sure C. Now touch the complex Hodge bundle over C with the Magic Wand and 
find its irreducible component containing the two integer cocycles responsible for 
the diffusion. Find or estimate the corresponding Lyapunov exponents and the 
problem is solved. This strategy evokes three problems: how to find an orbit clo¬ 
sure? How to find irreducible components of the Hodge bundle? How to estimate 
the top Lyapunov exponent of a given invariant flat subbundle? 

B.l. Classification of invariant suborbifolds. 

Problem 1. Classify all GL(2, M )-invariant suborbifolds in , d n ). 

The invariant suborbifolds in the strata of quadratic differentials are not men¬ 
tioned since replacing any flat surface corresponding to a quadratic differential 
by its canonical ramified double cover on which the induced quadratic differential 
becomes a global square of a holomorphic 1-form we get an associated invariant 
suborbifold in the corresponding stratum of Abelian differentials. 

In genus 2 the problem is solved by C. McMullen IMcMli . There are very 
serious advances in this problem for the stratum ??(4), due to Aulicino, Nguyen and 
Wright jANWI for H odd { 4) and Nguyen and Wright |ANW| for H hyp (4). Certain 
finiteness results for the number of primitive Teichmiiller curves are obtained by 
Bainbirdge and Moller jBMj and by Matheus and Wright |MaWrj . Situation in the 
Prym loci is very well understood due to the work of Lanneau and Nguyen ILN1I , 
lLN2l and lLN3l . 

The Theorem of Eskin-Mirzakhani-Mohammadi states that such suborbifolds 
correspond to complex affine subspaces in cohomological coordiantes. Due to result 
of Avila-Eskin-Moller |AEMo| the projection of such affine subspace to absolute 
cohomology should be a symplectic subspace: the restriction of the intersection form 
to this subspace is nondegenerate. Results of A. Wright IW 11 impose conditions 
on the field of definition. At some point all these constraints started to seem so 
restrictive that there was a doubt whether any really new invariant suborbifolds 
are left? Namely, there is a separate question of classification of Teichmiiller curves 
(both arithmetic and non arithmetic), and there are plenty invariant suborbifolds 
which can be obtained from Teichmiiller curves or from connected components of the 
strata by various ramified covering constructions. The question was whether there 
is anything else (with exception for some sporadic series of invariant suborbifolds 
specific for very small genera). 

A recent example of an invariant submanifold of enigmatic origin in 1-1(6) (which 
does not fit any of the known schemes mentioned above) was recently found by 
M. Mirzakhani and A. Wright. For a brief overview of the state of the art in this 
direction see |W3I . 

B.2. Classification of invariant subbundles of the complex Hodge bundle. 

The complex Hodge bundle H over the moduli space M g of curves has the ho¬ 
mology space iL 1 (C;C) as a fiber over the point represented by the curve C. This 
bundle can be pulled back to any stratum ‘ H(m \,..., m n ). 

Problem 2. Find the decomposition of the complex Hodge bundle H ^ over any 
given GL(2, M.)-invariant suborbifold into irreducible GL(2,R )-equivariant subbun¬ 
dles. 
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The fact that such decomposition exists is proved by S. Filip jFilj . Here “irre¬ 
ducible” should be understood in the sense that there is no further splitting even 
when we pass to any finite ramified cover of the GL(2, R)-invariant suborbifold. 

The problem is meaningful already for the strata! It is absolutely frustrating, 
but we do not have a proof that the only equivariant subbundles of H ^ over a con¬ 
nected component of any stratum is the tautological subbundle and its symplectic 
complement. We do not know either whether H\ over any connected component 
of any stratum of quadratic differentials is irreducible in this sense. 

The current tools allow, in principal, to prove the latter two facts more or less by 
hands for some low-dimensional strata. Namely, one can start with an arithmetic 
Teichmiiller curve, compute certain number of monodromy matrices and then use 
technique of |MMoYj or the method of Eskin as in |FMZ| to prove irreducibility of 
the complex Hodge bundle over the Teichmiiller curve. As a consequence of |Fil| 
one gets the irreducibility over over the ambient stratum. However, what is really 
needed is some general proof for all strata at once. 

A related question is 

Problem 3. What groups are realizable as Zariski closures of leafwise monodromy 
groups of equivariant irreducible blocks of the complex Hodge bundle H ^ over GL(2, R) - 
invariant suborbifolds ? 

The original guess of Forni-Matheus-Zorich [FMZI states that this group is al¬ 
ways SU (p, q ) for appropriate p and q. The paper of S. Filip |Fi2| shows that general 
Hodge-theoretical arguments admit a priori larger list (including some more sophis¬ 
ticated representations of SU (p, q)). However, it is not clear which of these groups 
(representations) are realizable as Zariski closures of leafwise monodromy groups 
of equivariant subbundles of the complex Hodge bundle over GL(2 ,WL)-invariant 
suborbifolds (and not just over some flat subbundles of the complex Hodge bundle 
over abstract submanifolds of the moduli space). 

B.3. Estimates for individual Lyapunov exponents. Paper |EKZ2| provides 
a formula for the sum of the positive Lyapunov exponents of the complex Hodge 
bundle over along the Teichmiiller geodesic flow. Though there is no reason to hope 
for exact values of individual Lyapunov exponents, paper |Yul of Fei Yu conjectures 
that partial sums of Lyapunov exponents might be estimated through Chern classes 
of holomorphic vector bundles over Teichmiiller curves normalized by the Euler 
characteristics of Teichmiiller curves. More generally (though less precisely): 

Problem 4. Study extremal properties of the “curvature” of the Lyapunov sub¬ 
bundles compared to holomorphic subbundles of the Hodge bundle. Estimate the 
individual Lyapunov exponents. 

For example, estimates for A) 1- over hyperelliptic locus in the principal stratum 
and estimates for Af over the entire principal stratum Q{ 1,..., 1) of holomorphic 
quadratic differentials would provide estimates for the diffusion rate in certain fam¬ 
ilies of wind-tree billiards. 

B.4. Siegel—Veech constants in terms of an adequate intersection the¬ 
ory. The sum of positive Lyapunov exponents of the complex Hodge bundle over 
an SL(2,R)-invariant suborbifold C in a stratum ..., m n ) is expressed 
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in 1EKZ2] as 


Ai H- + X g — 


1 ^ “l” 2) 

12 ^ rrii + 1 

Z— 1 


* 


m-. 


where c area (C) is the Siegel-Veech constant of £. Currently there are two formulae 
for c area (C) for two extremal cases of C. When £ is a connected component of a 
stratum, the Siegel-Veech constant c areo (£) is expressed as a polynomial in vol¬ 
umes of simpler principal boundary strata (normalized by the volume of the initial 
stratum); see lEMsZi for the strata of Abelian differentials and |Gj| for the strata 
of quadratic differentials. When £ is a Teichmiillcr curve, c area (£) is expressed as 
the integral of the Chern class of the determinant bundle over £ normalized by the 
Euler characteristic of £, see |l\ 1 . [BwMol . j EKZl ). The challenge is to construct 
a bridge between these two cases: 


Problem 5. Express c area (£) in terms of an appropriate intersection theory. 


Here we do not mean some kind of asymptotic limit formulae, but something in 
the spirit of ELSV-formula for Hurwitz numbers, see lELSVI . 

There is certain resemblance between the “hyperbolic regime” studied in [Mill 
|Mi3j by M. Mirzakhani and the “flat regime” studied in IEMsZ] . M. Mirzakhani 
used dynamics on moduli space to relate the length functions of simple geodesics 
on hyperbolic surfaces to the Weil Peterson volumes of the moduli spaces A4 9jn of 
punctured Riemann surfaces, and also to relate the Weil Peterson volumes to the 
intersection numbers of tautological line bundles over Mg t n- 

Morally, we have somehow similar situation in a parallel flat world. The step of 
relating the counting functions for simple flat geodesics (for the flat metrics in the 
same conformal class as the original hyperbolic metric) to polynomial in volumes of 
the strata with respect to Masur-Veech volume form is already performed in IEMsZ] 
and m■ The challenge is to accomplish the second step and to relate these volumes 
to an adequate intersection theory. 

Actually, certain parallel between hyperbolic and flat word manifests in further 
aspects. For example, there are conjectural simple asymptotic formulae for large 
genera for the Weil Peterson volumes |Mi3j and for Masur-Veech volumes. About 
ten years ago A. Eskin and one of the authors conjectured a very simple and ex¬ 
plicit asymptotic formula for the Masur-Veech volume. For the principal stratum 
H.( l 2 ® -2 ) this conjecture was recently proved in [CMoZj by D. Chen, M. Moller, 
and D. Zagier. 


B.5. Dynamics on other families of complex varieties. One more challenging 
direction of study is dynamics of the complex Hodge bundle over geodesic flows over 
moduli spaces of higher-dimensional complex manifolds. 

Problem 6. Study dynamics of the Hodge bundle over geodesic flows on other 
families of compact varieties. Are there other dynamical systems (compared to 
billiards in rational polygons) which admit renormalization leading to dynamics on 
families of complex varieties ? 

Some experimental results for families of Calabi Yau varieties are recently ob¬ 
tained by M. Kontsevich |K2| . S. Filip studied in |Fi3j families of K3-surfaces. 
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